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A broadband and broad-angle low- scattering metasurface is designed, fabricated, and characterized. Based 
on the optimization algorithm and far-field scattering pattern analysis, we propose a rapid and efficient 
method to design metasurfaces, which avoids the large amount of time-consuming electromagnetic 
simulations. Full-wave simulation and measurement results show that the proposed metasurface is 
insensitive to the polarization of incident waves, and presents good scattering-reduction properties for 
oblique incident waves. 



Metamaterials are engineered materials with unusual constitutive parameters such as negative permittivity 
and permeability 1 " 4 . They are usually composed by periodic resonant or non- resonant units with the 
characteristic dimensions much smaller than the wavelength 5 " 7 . One important application of the 
metamaterials is to achieve electromagnetic invisibility or transparency, which can greatly reduce the radar cross 
section (RCS) of metallic or dielectric targets 8 " 11 . The approach of optical transformation provides an efficient way 
to bend the incident waves along the boundary of a region to be hidden without perturbing the exterior fields 12 " 14 . 
Therefore perfect invisibility can be achieved with no backward or forward RCS. However, the objects inside are 
completely isolated, which makes such devices hindered in practical applications. Another promising method is 
based on the scattering cancellation 15 " 19 , where metasurface has been utilized to wrap the objects, and it can 
provide additional surface reactance in order to generate anti-phase scattered fields. But the bandwidth of the 
metasurface is usually very limited. The radar absorbing material can also be used for RCS reduction, which 
transforms electromagnetic energy into heat. But the change of the object temperature will inevitably increase the 
detection possibility by infrared detectors. 

In this paper, we develop a low RCS metasurface using the windmill-shaped unit. The phase of each unit is 
chosen to be distributed randomly. The reflected waves from each part of the metasurface will experience 
destructive inference due to the random phase difference, leading to electromagnetic diffusion for the scattering 
waves. The design of the proposed metasurface seems to be the counterpart of the traditional reflectarray, which 
requires the reflected waves to be steered over a range of directions. In order to redistribute the scattering energy 
into all the directions, a large phase swing (i.e., the available phase shift of the composing unit) is usually required, 
while the magnitude of the normalized reflection efficient for the unit is kept to approach zero (dB) at the same 
time. Similar techniques for RCS conduction based on the objects' coating have been reported in Refs. [15-17], 
where both the perfect conductor and artificial magnetic conductor (AMC) have been used to create a chess- 
board-like structure. The reflected waves from the two units have 180.phase difference, hence the destructive 
interference between the two waves will dramatically reduce the scattering field level. The bandwidth of such 
coating is also restricted since it is usually hard to design broadband AMC structures. The mechanism of the 
proposed metasurface mainly depends on the random distribution of the reflection phases, which is totally 
different from the mechanisms mentioned above. 



Results 

Metasurface is a two-dimensional equivalent of three -dimension metamaterial (see Fig. 1), which could be used to 
control the reflection and transmission properties for the incident waves 20 " 22 , and therefore achieve some inter- 
esting functionalities such as focusing and wave bending 23 " 25 . Usually a metasurface is composed by periodic or 
aperiodic elements, which are printed on the surface of a thin substrate. 
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Figure 1 | Reflection and transmission of electromagnetic waves on the 
metasurface. 

For a metasurface composed by large number of units with the 
same type and different geometrical dimensions, electromagnetic dif- 
fusions will be caused, resulting in low backward RCS 26 . However, it is 
usually hard to simulate the far-field scattering radiation pattern due 
to the burden of computation time and memory sources. Therefore 
the optimization of such metasurface for RCS reduction becomes very 
inefficient and impractical based on the simulation results. Actually 
the scattering pattern of such complex array can be rapidly synthe- 
sized when the mutual coupling between the adjacent elements are 
achieved through simple numerical simulations 27 " 28 . Let us consider a 
two dimensional metasurface composed by M*N units. When illu- 
minated by a plane wave, the total scattering wave from the metasur- 
face can be regarded as the superposition of the scattering wave from 
each basic element, which can be expressed as follows, 



Etotal= ^ y^Em,n(6,<p)' 



(i) 



where E m>n is the vectorial far-field scattered by the element [m,n]. 6 
and cp are the polar angle, and azimuthal angle, respectively. To obtain 
the scattering field of each element, the infinite periodic array model is 
adapted to calculate the coefficient E m>n in Eq. (1). We assume that 
every element is surrounded by the same cells when considering the 
influence of mutual coupling. This is an approximate model which 
does not correspond to the realistic metasurface configurations. 

Here the induced surface currents from each element with peri- 
odic boundaries can be extracted from the numerical simulations. 
Take a small metasurface of 5*5 units for example, (as shown in 
Fig. 1(a)), the metallic windmill-shaped structure is selected as the 
basic element, and a perfect conductor layer is attached to the back of 
the substrate. From the equivalence principle, there are three kinds 
equivalent current on the surface of each element: the equivalent 
conducting current J c on the conductor area and the equivalent elec- 
tric current J d and magnetic current M d on the dielectric area 29 " 30 . In 
our simulation, the full-wave electromagnetic simulation tool FEKO 
is used to extract the equivalent currents mentioned above. When the 
metasurface is illuminated by a plane wave, the scattering field of 
each element can be written as 



E s (x) = L c [J c ]+L d [h]+L d [M d 



(2) 



where 

L c [J c } = -j<j)A c - 
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dielectric 



and G(r : r') is the green function in free space. 

Here the dimension of example is 60 mm X 60 mm, with several 
random-sized windmill -shaped structures. The radiation pattern the 
metasurface in Fig. 2(a) is also simulated by FEKO to check the 
agreement between the synthesis method and the numerical method, 
as shown by the black line and the red line in Fig. 2(b). It is clear that 
the synthesis method in Eq. (1) provides satisfactory prediction of the 
scattering pattern of the metasurface from the comparison. There are 
some discrepancies in the sidelobes from the two results in Fig. 2(b), 
which is mainly due to the mismatch of the boundary conditions 
used in the simulations. From the comparison between the simulated 
and synthesized radiation pattern, both results agree well which 
means that this model sufficiently provides a powerful means for 
radiation pattern prediction and optimization. 

To get best performance of the metasurface, a particle swarm 
optimization (PSO) algorithm is utilized together with far-field pat- 
tern predition to achieve the optimal arrangement of the meta- 
atoms. In PSO algorithm every particle can be seen as an individual 
who searches for the optimum solution through cooperation in indi- 
viduals and sharing information, which is an intelligent optimization 
algorithm proposed by Kenndey and Eberhart 31 . It has been widely 
applied to the optimization of various electromagnetic problems due 
to the merits like simple principle, few parameters and fast conver- 
gence. The potential solution can be deemed to a particle, whose 
mathematical description of the speed and location of next genera- 
tion could be expressed as, 



J+ 1 



^4+ C\ x randi x (^pbest^ ik — + c 2 x rand 2 x (gbes^ ■ 



x ik ' 



-4) 

(9) 
(10) 



Where co is the inertia coefficient, c x and c 2 are the acceleration 
constants, which are used to adjust individual speed. rand Y and rand 2 
are random numbers between 0 and 1. pbes^ is the individual 
optimal location, gbest j is the global optimal location. In the design 
of metasurface, x denotes the size of basic unit, and the performance 
of metasurface is evaluated by fitness(x), which is described as, 

fitness = c 3 x max(£ 5 ) + c 4 x £ 5; x (11) 

where c 3 and c 4 refer to weight coefficients. E s denotes the total 
scattered field, E S} i denotes scattered field in the specific direction. 
Such fitness function aims at reducing the scattering energy mostly in 
the direction of maximum radiation. 
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Figure 2 | (a) A metasurface containing 5*5 elements, (b) RCS of the metasurface in Figure 2(a) from the numerical simulation and far-field pattern 
analysis. 



The design flow of the low RCS metasurface is shown in Fig. 3, 
where two modules have been used to get the desired optimal design. 
The PSO module evaluates the fitness, updates the particle speed and 
the population location (i.e. the meta-atom arrangements) in each 
iteration, and then sends the information to the far-field scattering 
pattern module. The latter gives the scattering pattern of the opti- 
mized metasurface based on the equivalent currents mentioned in 
last sections and calculates the reduction of RCS in the main scatter- 
ing directions. Then the value of fitness is computed and returned to 
PSO module. After some iterations, we can get the optimal arrange- 
ment of metasurface with lowest RCS as required. The kernel of this 
flowchart is that the PSO module determines various combinations 
of basic unit, whose performances are judged by the far-field pattern 
analysis module to find the best solution. 



A low RCS metasurface is designed, as shown in Fig. 4(a). The 
windmill shaped structure is selected to be the basic element. It is 
composed by an inner square ring and an outer irregular ring, and 
their resonance frequencies are designed to be very close in order to 
increase the operation bandwidth and the reflection phasing range. 
The curves of the reflection phase from 6 GHz to 14 GHz is almost 
linear, and they are nearly parallel with the change of L (the length of 
the branch of the outer ring), showing excellent broadband fre- 
quency response as required 32 . Here we choose the dimensions of 
the basic element in Fig. 4(a) asa=12 mm,g=0.4 mm, <i=0.55 mm, 
h=4 mm. The substrate used in the simulation and experiment is 
F4B (s r =2.65+i.001) with the thickness 4 mm. The reflection phase 
of the basic element is extracted through the numerical simulations 
with the change of L at 8 GHz, 10 GHz, 12 GHz and 14 GHz (see 
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Figure 3 | Flowchart for the design of low RCS metasurface. 
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Figure 4 | (a) Basic element of the low RCS metasurface and the layout of the designed metasurface, (b) the reflection phase from 6 GHz to 14 GHz with 
the change of L in Fig. 4(a), (c) the reflection phase at different frequency with the change of L. 



Fig. 4(c)). These curves are nearly parallel with each other, showing 
excellent linear frequency response, which is critical to guarantee the 
working bandwidth. It is clear that a large phasing range can be 
achieved when tuning the branch of the outer ring, which is import- 
ant to ensure the performance of the low RCS metasurface. 

In this design, for the simplicity of optimization, we choose ele- 
ments with discrete phase distributions instead of continuous ones as 
mentioned above. Here eight windmill- shaped elements with differ- 
ent dimensions have been selected as the composing unit of the 
metasurface. When the length of the branch L is 0.6 mm, 



0.825 mm, 1.05 mm, 1.26 mm, 1.47 mm, 1.75 mm, 2.05 mm and 
2.38 mm respectively, the reflected phase is 45 *n (n= 1 to 8) degrees 
at the central frequency f = 10 GHz, as illustrated in Fig. 4(b). 

The whole metasurface is made up of 144 windmill-shaped units, 
and we need to find the optimal combination for the desired low RCS 
metasurface based on the eight basic elements. In our algorithm, a 
population size of 400 and a maximum of 100 iterations are prede- 
fined. It will take 9 hours to finish the every iteration based on the 
numerical simulation. However, the time spent in every iteration 
using the far-field scattering module is approximately 10 minutes 
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Figure 5 | (a) Far-field scattering pattern (<p = 0°) from numerical simulation (by FEKO) and far-field scattering pattern analysis, (b) simulated RCS 
reduction for normal incidence from 6 GHz to 14 GHz. 
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Figure 6 | (a) The experiment setup for RCS reduction, (b), (c) Simulated and experimental RCS reduction for various incident angles with vertical 
polarization, (d), (e) Simulated and experimental RCS reduction for various incident angles with horizontal polarization. 



with the same hardware, which is much more efficient than before. 
After some iterations, the arrangement of the basic elements is finally 
determined with the layout shown in Fig. 4(a). 

Fig. 5(a) shows the scattering pattern of the metasurface from full- 
wave simulations and far-field pattern prediction analysis. Both 
results seem to agree very well, which validates the whole design flow 
in Sec. 3. The RCS reduction for a bare metallic plate and the designed 
metasurface at different frequencies is illustrated in Fig. 5(b), where 
the RCS reduction is greater than 10 dB from 7 GHz to 13 GHz. 

Discussion 

In order to verify the theoretical analysis above, a larger metasurface 
containing 16 by 16 units, optimized by the algorithm, has been 



fabricated and tested, as depicted in Fig. 6. The transmitting and 
receiving horn antennas have been connected to a vector network 
analyzer (Agilent N5230C) to measure the RCS of the metasurface. 

The simulated and measured RCS reductions for vertical and 
horizontal polarizations are shown in Figs. 6(c) (d) and Figs. 6(e) (f) 
respectively. The maximal RCS reduction can be found to exceed 
30 dB as we can see from the measured results. The metasurface 
shows excellent broadband properties, where the RCS reduction is 
over 10 dB from 7 GHz to 14 GHz. With the increase of the incident 
angle, the bandwidth decreases a little due to the phase aberrations. 
Similarly, the maximum of RCS reduction is also decreased. 
However, low RCS properties are still kept for the proposed 
metasurface. 
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We have proposed the design algorithm for a low RCS metasurface 
with ultra-low backward scattering. The optimization method is used 
together with the far- field scattering pattern technique to find the 
desired broadband metasurface. The simulated and measured results 
demonstrate that the metasurface can effectively decrease the back- 
ward scattering from 7 GHz to 14 GHz, which is important for 
stealth applications in the future. 

Methods 

The sample is placed as high as the horn antennas in the experiment. The distance R 
(3 m) between antennas and sample is chosen far enough to avoid the near field effect. 
Pyramidal absorbing materials have been placed around the sample to decrease the 
unwanted reflections from the surroundings. In addition, a copper plate with the 
same size has also been measured for comparison. In the case of normal incidence, the 
transmitting and receiving horn antennas are placed adjacently, where the real 
incident angle is about 3 degree because of the restricted size of horn antennas. Both 
the transmitting and receiving horn antennas can move along the circumference trace 
to obtain RCS reductions at different scattering angle 0 scatt , where we only consider 
the case that the scattering angle 0 scatt equals to the incident angle 0 inc . 
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